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ABSTRACT 

Solutions  have  been  obtained  for  the  temperature  and  moisture  distribu- 
tions within  the  atmospheric  boundary  layer  of  an  axisymmetric  hurricane 
model.   The  intensity  of  the  hurricane  is  related  to  the  equivalent  potential 
temperature  gradient  produced  by  a  balance  of  heat  sources  and  advection 
within  the  boundary  layer.   Solutions  are  obtained  using  the  bulk  aerodynamic 
transport  equations  or  applying  a  two- layer,  baroclinic  boundary  layer  model 
by  Cardone  (1969).   Equilibrium  maximum  wind  speeds  vary  from  minimal  hurri- 
cane force  to  about  80  m/sec  for  fixed  ocean  temperatures  between  27C  and  31C. 
Observations  of  the  oceanic  heat  loss  in  actual  storms  are  necessary  to  estab- 
lish the  ratio  of  heat  flux  to  momentum  transfer  in  the  high  wind  speeds 
treated  by  the  model.   Various  applications  of  the  model  are  proposed,  with 
an  ultimate  goal  of  a  time-dependent  simulation  of  hurricane-ocean  coupling. 

INTRODUCTION 

A  number  of  axisymmetric  vortex  models  (e.g.,  Yamasaki,  1968a  and  1968b; 
Ooyama,  1969;  Rosenthal,  1971)  have  been  developed  to  describe  the  evolution 
of  the  intense  wind  region  characteristic  of  hurricanes.   The  specified 
initial  conditions  are  mainly  designed  to  insure  low-level  inflow,  so  that 
the  vortex  will  eventually  become  unstable,  and  are  not  related  to  observa- 
tions.  As  a  result,  a  comparatively  long  "organizational  period"  is  required 
before  deepening  of  the  vortex  ensues.  A  primary  result  of  these  models  is 
the  very  realistic  simulation  of  the  wind  structure  near  the  center  of  the 
hurricane.  Because  the  initial  conditions  are  arbitrarily  specified,  rather 
than  from  real  data,  the  intermediate  stages  of  the  vortex  simulation  may  not 
be  realistic.  These  models  do  not  allow  forcing  by  external  circulations. 


Instead,  the  philosophy  is  to  displace  the  exterior  boundaries  of  the  model 
as  far  as  necessary  to  assure  that  the  boundary  conditions  will  not  affect 
the  solution  near  the  center.   A  complete  assessment  of  external  circula- 
tions, and  other  asymmetric  effects,  on  hurricane  simulation  must  await  fur- 
ther development  (see  Anthes  et  al,  1971)  of  three-dimensional  models  started 
with  real  initial  data. 

The  objective  of  this  paper  is  to  describe  a  simplified  model  of  the 
properties  within  the  hurricane  boundary  layer.  Although  the  ultimate  pur- 
pose is  to  simulate  the  time  variation  of  these  properties  due  to  changes  in 
sea-surface  temperatures,  in  the  present  paper  the  emphasis  will  be  on  the 
adjustment  to  different  fixed  sea-surface  temperatures.   It  is  well  known 
that  tropical  storms  develop  into  hurricanes  only  over  warm  tropical  oceans. 
Indeed  the  oceanic  surface  layers  must  provide  the  heat  necessary  to  produce 
the  unique  warm  core  of  the  mature  hurricane.  The  preferential  release  of 
the  heat  in  deep  convective  towers  near  the  hurricane  eye  results  in  central 
temperatures  much  higher  than  the  environment.  To  a  first  approximation, 
the  heat  gain  within  the  boundary  layer  may  be  linearly  related  to  the 

tangential  wind  speed  at  the  top  of  the  layer  (Vq),  but  the  frictional  loss 

2 
of  momentum  to  the  ocean  increases  as  Vq  .   Thus  a  given  radial  profile  of 

wind  speed  essentially  specifies  the  heat  gain  which  sustains  the  humricane, 
as  well  as  the  frictional  loss  which  tends  to  offset  further  increases  in 
wind  speed.   On  this  basis,  Riehl  (1963)  developed  relations  for  the  unique 
wind  profile  (and  associated  thermal  structure)  for  which  these  offsetting 
tendencies  are  in  balance;  that  is,  for  the  conditions  in  which  a  steady 
state  hurricane  should  result.  However,  an  actual  storm  may  not  reach  its 
full  potential,  or  may  be  in  steady  state  for  only  a  short  time,  because  the 
characteristics  of  the  ocean  and  surrounding  environment  vary  in  space  and 
time.   A  prime  motivation  of  this  study  was  to  develop  a  model  capable  of 
simulating  the  symmetric  component  of  actual  storms. 

2 


DESCRIPTION  OF  THE  MODEL 

The  empiricism  of  the  model  is  based  on  three  assumptions.   First,  the 

wind  profile  at  the  top  of  the  boundary  layer  is  specified  by  the  relation 

1/2  •  :  •■  ■' 

-    ■  Vgr    =  constant  ,  (1) 

outside  the  radius  of  maximum  wind  speed  (r.).   Between  the  vortex  center  and 
r.  the  winds  are  assumed  to  obey  VqT       =   constant,  which  represents  solid- 
body  rotation.   Justification  of  (1)  was  given  by  Gray  and  Shea  (1973)  who 
found  an  exponent  of  0.47  +  0.30  in  Vq  profiles  from  hurricanes  in  various 
stages  of  development.   Rather  than  using  individual  Vq  profiles,  Riehl  (1963) 
averaged  profiles  around  the  storm  and  found  the  exponent  (1/2)  applied  in 
many  cases,  particularly  in  quasi-steady  state  hurricanes.  Some  of  the  addi- 
tional variation  which  Gray  and  Shea  found  must  be  due  to  transient  features 
in  the  wind  field  as  the  storm  adjusts  to  varying  internal  or  external  condi- 
tions.  In  simulating  the  free-atmosphere  as  a  series  of  "equilibritim  states" 
represented  by  (1),  these  transient  features  are  assumed  to  be  of  less 
importance  than  the  symmetrical  circulation.   It  should  be  emphasized  that 
the  coefficient  in  (1)  was  determined  from  aircraft  data  only  within  about 
150  km  of  the  center.   This  is  the  region  in  which  the  symmetry  assumption  is 
most  justifiable,  and,  fortunately,  is  the  primary  region  of  interest  in 
evaluating  the  effects  of  the  hurricane.  Lack  of  systematic  observations  at 
larger  radii  prevent  an  accurate  description  of  the  outer  region  in  which 
the  external  effects  are  likely  to  be  manifested  in  actual  hurricanes.  There- 
fore, beyond  r   [defined  in  (5)]  the  exponent  in  (1)  was  adjusted  to  always 
reduce  the  wind  speed  at  2r  to  5  m  sec   .   Specifying  a  different  exponent 
in  the  outer  region  had  an  insignificant  effect  on  the  properties  near  the 
center  compared  to  experiments  in  which  (1)  was  applied  throughout  the  region. 


The  second  empirical  formula  relates  the  maximum  wind  speed  to  the  in- 
crease in  the  equivalent  potential  temperature  of  the  boundary  layer  (60  ) 
due  to  the  oceanic  heat  source.   Such  a  relation  was  proposed  by  Riehl  (1963) 
through  the  use  of  the  cyclostrophic  assumption  and  observations  relating 

the  surface  pressure  at  the  radius  of  maximum  wind  to  60  , 
^  e 

6p  =  -  2.5  60   ,  (2) 

where  the  6  symbols  indicate  deviations  from  p  =  1005  and  0  =  350K.   In  the 
present  model,  a  gradient  wind  expression 


=Vn(f+— )  (3) 


p  Sr    t)  ~    r 


was  used  instead  of  the  cyclostrophic  assumption.   Although  Ooyama  (1969) 

showed  that  the  inertial  term  is  important  near  the  surface,  this  term  may  be 

dropped  if  (3)  is  averaged  over  the  depth  of  the  boundary.   Upon  integration 

of  (3)  between  the  radii  over  which  (2)  applies,  the  tangential  component 

(vq  )  at  the  radius  of  maximum  wind  (r.)  is 
0.  1 
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0.5  r                  I             __  r  -r, 

Vn     =  -['2f(Vr.Y  -r.)+t/4f2(v/r.r     -  r.)2  +io6  _2 — i(e     .9 

0.        r     -r.                   X  o      1      If         ^10  1                    pr          e.      e 

101                                     »  010 


)](4) 


The  dominant  term  in  this  expression  is  that  involving  the  difference  in  0 

at  r.  and  at  some  radius  r  where  p  =  1005  and  0  =  350K,  which  follows 
1  o       *^  e       ' 

from  (2).  Bell  and  Tsui  (1973)  found  a  similar  relation  between  pressure 
and  0  changes  for  tropical  cyclones  near  the  Asian  coast,  but  with  a  co- 
efficient of  2.25.  Because  of  the  difficulty  in  establishing  the  desired 
surface  pressure  at  r.,  and  only  a  10^  percent  difference,  the  coefficient 
in  (2)  was  adopted.  It  is  clear  from  (4)  that  the  radius  of  maximum  wind 
(r . )  should  coincide  with  the  maximum  0  within  the  boundary  layer.   In  the 


present  model,  r,  was  chosen  as  the  radius  at  which  6  ,  averaged  over  five 
grid  increments  (Ar  =  3km),  was  a  maximum.   This  assumption  makes  r.  coin- 
cident with  the  eyewall  clouds,  because  the  boundary  layer  air  parcels  ' 
with  the  largest  heat  and  moisture  content  would  be  expected  to  be  most 
favorable  for  deep  convection.   Shea  and  Gray  (1973)  show  that  r,  tends  to 

be  coincident  with  the  eyewall  in  the  most  intense  hurricanes,  but  r   is 

i 

larger  than  the  eyewall  radius  in  weaker  storms.   The  radius  (r  )  at  which 

o 

p  =  1005  and  9  =  350K  may  be  taken  to  represent  an  external  influence,  and 
must  be  defined  in  terms  of  other  parameters.   As  the  p  and  6  gradients  are 
small  outside  the  central  region,  the  inaccuracy  in  specifying  r   is  a 
second-order  effect  compared  with  r..   For  convenience  we  choose 

r^  =  (f  ^e/i)  (5) 

1  ,  ^  ■ 

which  is  approximately  the  radius  at  which  Vn  tends  to  zero  in  the  outflow 
layer  of  the  storm,  if  the  air  originating  in  the  maximum  wind  region  con- 
serves absolute  angular  momentum.   Since  r   is  about  200km  for  f^(10   )sec 

and  typical  values  of  Vq  and  r.,  this  again  represents  the  region  of  appli- 

i 
cation  of  the  model.   A  schematic  representation  of  the  wind  structure  in 

the  atmospheric  vortex  is  given  in  Figure  1. 

The  above  relations  represent  an  empirical  model  of  the  wind  speeds  in 

hurricanes  as  a  function  of  the  6  distribution  within  the  boundary  layer. 

e 

Rather  than  determine  the  maximum  0   from  observations  as  proposed  by  Riehl 

e 

(1963),  the  0  fields  will  be  calculated  from  time- dependent  equations  for 
the  potential  temperature  (0)  and  specific  humidity  (q),  ^- 


0^  =   0   exp    (^) 
P 

^t  = 

S(rv   0)      A,    QN        gQ 
^     r          d(cD0)        ^^s 

r<^r       "^     Sp            C   Ap 
P 

,  and  (6) 


N      ^(rv  q)   ^.   .    gQ 

In  these  equations  C  and  L  are  the  specific  heat  and  latent  heat  of  vapori- 
zation.  Only  the  sjnmmetrical  circulation  is  considered,  so  that  the  radial 
(v  )  and  vertical  (cu)  components  of  motion  are  related  by 


-X      ^rv 
(3cD  _  _   r 

§p     r^r 


(8) 


with  the  further  approximation  that  oj  =  0  at  the  surface.  Radiative  exchanges 
are  assumed  to  be  negligible  and  the  sensible  (Q  )  and  latent  (Q  )  heat 
fluxes  are  assumed  to  decrease  to  zero  within  Ap  =  100  mb  of  the  surface. 
The  actual  vertical  heat  fluxes  through  the  top  of  the  boundary  layer  are 
determined  by  the  mean  flow  represented  by  03  and  convective  towers.   Normally 
the  convective  heat  flux  is  parameterized  in  terms  of  the  mean  flow;  however, 
this  resulted  in  supersaturation  in  the  present  model,  as  was  the  experience 
of  Rosenthal  (1971).   The  third  important  empirical  aspect  of  the  model  was 
the  formulation  of  a  critical  relative  humidity  profile,  such  that  the 
moisture  in  excess  of  the  critical  value  was  arbitrarily  removed  from  the 
boundary  layer.   The  convective  heat  flux  6o6  was  neglected. 

Because  the  radial  gradient  of  6   is  essentially  determined  by  the 
moisture  gradient,  the  specification  of  the  critical  relative  humidity  profile 
was  a  primary  factor  in  determining  the  intensity  of  the  model  hurricane. 
As  a  first  attempt,  the  critical  relative  humidity  for  the  air  in  the  exterior 
region  was  specified  as  the  value  which  would  produce  a  lifting  condensation 
level  (LCL)  of  60  mb  above  the  surface.   A  decrease  in  the  typical  cloud  base 
height  toward  the  center  is  commonly  observed  in  hurricanes,  although  the 
magnitude  of  the  decrease  has  not  been  documented  as  a  function  of  maximum 
wind  speed.   In  this  case,  a  linear  decrease  to  a  LCL  of  30  mb  at  the  eyewall 


was  assumed.   Profiles  other  than  a  linear  decrease  were  also  used,  but 
these  did  not  change  the  results.  These  LCL  values  corresponded  to  external 
and  internal  critical  relative  humidity  values  of  about  77  and  88^,  respec- 
tively, for  the  range  of  sea  surface  temperatures  used  here.  These  values 
seem  realistic  as  one  might  expect  the  external  air  to  be  somewhat  drier 
because  of  subsidence,  and  the  interior  air  to  be  more  moist  than  normal 
(about  80^  in  the  subcloud  layer).   In  the  second  approach,  the  critical 
relative  humidity  values  were  arbitrarily  specified  as  a  linear  profile  be- 
tween 95^  inside  r.  and  75^^  outside  r  .   These  values  might  be  taken  as 
representative  of  the  maximum  conditions  one  might  expect,  and  the  results 
for  this  specification  will  be  referred  to  as  the  maximum  hurricane.  A 
95/^  value  implies  a  cloud  base  of  about  150  m  near  the  center. 

BOUNDARY  LAYER  MODELS 

The  surface  fluxes  Q  ,  Q  and  T^  required  to  close  the  solutions  are 

s   e      6s   ^ 

determined  by  the  boundary  layer  model.   For  fixed  oceanic  conditions,  the 
model  will  tend  to  steady  state  as  the  radial  and  vertical  advection  of 
heat  and  moisture  within  the  boundary  layer  balance  the  oceanic  sources 
Q  and  Q  ,  respectively.   In  this  formulation  the  air-sea  temperature  and 
humidity  differences  will  vary  with  radius  as  determined  from  (6)  and  (7). 
The  first  approach  to  the  boundary  layer  was  to  assume  the  bulk  aerodynamic 
equations  for  momentum,  heat  and  moisture  fluxes: 


Cp  =  I  C„  =  1.1  X  10"^  +  4  X  10"^  |vj 


where  C„  is  the  heat /moisture  exchange  coefficient,  Tq  is  the  surface  shear- 

H  OS 

ing  stress,  q   is  the  specific  humidity  corresponding  to  the  saturation  vapor 
pressure  at  the  sea  surface  temperature  T  ,  and  the  total  wind  speed  (V  )  is 
0.7  times  the  wind  speed  at  the  top  of  the  boundary  layer.   If  the  shearing 
stress  vanishes  at  the  top  of  the  boundary  layer,  the  steady-state,  tangential 
equation  of  motion  may  be  used  to  determine  the  mean  radial  flow  within  the 
layer, 

V   =   ^   ;x .  (12) 

OVn     Vn 


r 


(f  +  v-^  +  — ) 

dr    r 


The  set  (1)  -  (12)  is  thus  closed  for  a  given  sea  surface  temperature,  and 
may  be  solved  as  an  initial  value  problem. 

As  an  alternative  to  the  bulk  aerodjmamic  equation  approach,  a  two-layer, 
baroclinic  model  developed  by  Cardone  (1969)  was  adopted.  The  Cardone  model 
of  the  marine  boundary  layer,  hereafter  referred  to  as  CBL,  was  an  extension 
of  the  neutral,  fixed-terrain  model  described  by  Blackadar  (1965).  The  two- 
layer  CBL  model  is  schematically  illustrated  in  Figure  2.   In  the  surface 
layer  the  eddy  viscosity  (K^)  is  a  function  of  height  and  atmospheric  sta- 
bility, according  to  the  Monin-Obukhov  similarity  theory.   Specifying  a 
constant  K^,  equal  to  its  value  at  the  top  of  the  surface  layer,  and  replac- 
ing f  by  (f  +  Vn/r),  results  in  a  modified  Ekman  spiral  in  the  upper  layer. 
Following  Cardone  (1969),  the  surface  roughness  (Z  )  in  meters  is  internally 
determined  as  a  function  of  the  friction  velocity  (u^^^), 

Z  =  (>'^^(10     )  ^  4.28(10"^)u^^  -  4.43(10"'^)        (13) 
o      u^ 

The  desired  boundary  layer  wind  profiles  and  associated  heat  flux  can 
be  determined  from  the  solution  of  the  following  set  of  simultaneous 
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equations. 


where 


C^  =V2  k  sin    (J  -  a^)/[ln  B^R^   -   ta^)]  (U) 

S^=     2  k  B^   sin^  CLJ   (t)(L^)  (15) 


'J 


Cp  =  '^*/V   >    the  geostrophic  drag  coefficient, 

k  =  0.4,   von  Karman's  constant, 

OC     =  surface  cross-isobaric  angle, 

R     =  V  /fZ    ,    surface  Rossby  number, 
o  g       o  ■'  * 

(|)(L^)   =  non-dimensional  wind  shear   function, 
i|f(L^)   =  wind  profile  function, 
L^  =  h/L',    stability  index, 

h  =  B^V   /[f  +  Vg/r] ,   height  of  constant   flux  layer, 

-4 
B     =  3  X  10      ,   dimensionless   constant, 
o  ' 

L'    =  a^"^  uj-    e[ln(10/Z^)-   i|;(10/L' )  ] /[gk^CT     -T      )  ] ,    the 
•  n  "  o  va     vs 

gradient   form  of  the  Monin-Obukhov   length  applied  at   10  m, 

a,    =  K^/Kw  ratio  of  the  heat  transfer  coefficient  to  the  eddy 

viscosity, 

9  =  mean  potential  temperature  in  the  boundary  layer,  and 

T  -T   =  virtual  temperature  difference  between  10  m  and  the  sea 
va  vs  ^ 

surface. 

Expressions  for  <P(Lj^)  and  tCL^)  and  the  iterative  method  for  solving  the 

highly  implicit  set  (13)  -  (15)  are  given  by  Cardone  (1969).   Solutions  for 

u^,  a  and  L'  specify  the  desired  wind  profiles  and  implied  heat  flux  in 

the  constant  flux  layer.   It  should  be  noted  that  the  virtual  temperature 

difference  (T  -T  )  is  used  to  represent  both  sensible  (Q  )  and  latent 
va  vs  s 

(Q  )  fluxes  from 


H  =  Qg  +  Q^  =  [-u^^  CpP  a^   S]/[kg  L']        (16) 


Individual  Q  and  Q   fluxes  are  specified  according  to  the  ratio  (T  -T) 
s      e  w 

/(^L  ~^)»   The  radial  wind  component  was  determined  by  integrating  the  wind 
profile  given  by  the  similarity  theory  at  four  points  within  the  constant 
flux  layer  and  by  the  modified  Ekman  spiral  above. 

Use  of  the  CBL  in  the  atmospheric  vortex  assumes  applicability  of  the 
similarity  theory  in  wind  regimes  for  which  observations  are  lacking.  Thus 
the  model  results  must  be  critically  examined;  hopefully,  the  results  will 
stimulate  interest  in  devising  and  obtaining  the  measurements  necessary  to 
improve  the  boundary  layer  representation  under  hurricane  conditions. 

Computational  aspects.   The  combined  vortex -boundary  layer  model  was 
solved  as  an  initial  value  problem  with  centered  space  and  time-differences 
after  an  initial  forward  step.  Forward  time  steps  were  also  used  periodi- 
cally to  reduce  the  time  separation  of  the  solutions.  Ocean  temperature 
values  were  constant  in  space  and  time.   Because  the  model  is  to  be  started 
from  non-zero  initial  wind  speeds,  temperature  and  moisture  fields  that  are 
consistent  with  the  vortex  intensity  must  be  specified.   Although  an  esti- 
mate of  9  along  the  radius  may  be  obtained  from  (2)  -  (4)  using  the  Vq 
e  0 

profile,  the  values  of  9  (or  T)  and  q  must  also  be  consistent  with  the  heat/ 
moisture  source  and  the  advection.   As  sufficient  data  were  not  available  to 
specify  all  of  these  variables,  two  techniques  were  used  to  specify  initial 
values.   In  the  first  experiments  the  relative  humidity  was  specified  and 
the  Clausius-Clapeyron  equation  and  the  9  definition  were  solved  iteratively 
for  the  temperature.  An  alternative  approach  was  to  relate  the  atmospheric 
temperature  to  the  oceanic  temperature,  and  to  determine  the  specific  humidity 
from  the  9  value  calculated  from  (2)  -  (4).   Thus  in  regions  where  the 
initial  air  temperature  was  too  high,  the  initial  q  field  was  too  low.   If 


10 


the  ocean  temperature  were  fixed,  the  radial  temperature  and  humidity  profiles 
approached  the  same  values  within  a  few  hours  with  both  initialization  tech- 
niques.  Consequently,  the  second  alternative  was  adopted  by  specifying  an  air 
temperature  IC  less  than  the  ocean  temperature. 

RESULTS 

Effect  of  Ocean  Temperature 

A  comparison  of  the  maximum  tangential  wind  speed  for  the  two  critical 
relative  humidity  specifications  is  shown  in  Table  1  for  the  two  boundary  layer 
versions.   In  both  cases  there  is  a  systematic  increase  in  maximum  Vq  as  the 
sea  surface  temperature  is  increased.   Ooyama  (1969)  also  found  a  change  of 
about  20  m  sec   for  a  2C  change  in  sea  temperature.  When  the  moisture  speci- 
fication dependent  on  the  cloud  base  is  used,  minimal  hurricane-force  winds  are 
obtained  for  an  ocean  temperature  of  28C.  With  the  maximum  relative  humidity 
specified  as  95^,  hurricane- force  winds  are  obtained  with  a  sea  temperature  of 
27c.   Similar  differences  in  maximum  wind  speed  are  obtained  at  the  same 
ocean  temperature  for  the  two  moisture  specifications,  irrespective  of  the 
boundary  layer  model.  The  differences  caused  by  the  moisture  specification 
decrease  slowly  with  increasing  sea  surface  temperature. 

Riehl  (1963)  noted  the  threshold  0  value  for  producing  surface  pressures 

e 

capable  of  sustaining  tropical  storms  and  hurricanes  was  350K.   Parcels  with 

surface  properties  corresponding  to  9  >  350K  are  required  for  deep  convection 

and  warming  of  the  tropospheric  column.   For  higher  oceanic  temperature, 

the  air  temperatures  and  specific  humidities  throughout  the  boundary  layer 

are  increased.   Thus  the  environmental  9  values  (at  r  )  tended  to  increase 

e  o 

along  with  the  values  near  the  eyewall.  An  increase  in  9  at  r.  relative  to 

9  at  r  is  associated  with  a  pressure  decrease  as  specified  in  (2),  and 
e     o 

also  with  the  maximum  wind  through  (4).   However,  the  9  of  the  environmental 
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air  exceeds  350K  for  the  experiments  with  high  ocean  temperatures.   If 

only  the  increase  in  9  over  the  value  at  r  is  assumed  to  contribute  to 
■^  e  o 

the  surface  pressure  drop,  the  maximum  wind  speeds  shown  in  parentheses 
in  Table  1  result.   Only  a  small  increase  in  wind  speed  is  predicted  for 
water  temperatures  above  30C.  Within  the  framework  of  three  basic  assump- 
tions of  the  model,  such  an  interpretation  would  suggest  a  plateau  of 
maximum  intensity  of  storms  regardless  of  water  temperature  in  excess  of 
a  certain  value.   For  a  95^  specification  of  the  critical  relative  humidity 
at  the  eyewall,  the  66  values  are  about  13-15K  over  the  entire  range  of 
ocean  temperatures  in  Table  1.  Wind  speeds  exceeding  the  plateau  value 
could  only  be  expected  if  anomalously  cool  and  dry  conditions  existed 
around  a  tropical  storm  over  a  warm  ocean.  Such  conditions  might  be  ex- 
pected as  the  storm  nears  a  coastal  region,  if  the  inflowing  boundary  layer 
air  originating  over  land  has  a  sufficiently  long  trajectory  over  water 
before  entering  the  central  regions  of  the  storm.  One  can  speculate  that 
the  combination  of  a  favorable  environment  and  the  warm  coastal  waters 
could  lead  to  intensification  near  the  coast,  prior  to  the  filling  phase 
after  the  storm  center  crosses  the  coast. 

It  is  somewhat  satisfying  that  the  crude  representation  in  terms  of 
the  bulk  aerodynamic  transport  equations  produces  maximum  intensities  which 
are  quite  similar  to  those  obtained  from  the  more  sophisticated  CBL.  The 
CBL,  incorporating  the  empirical  constants  as  specified  by  Cardone  (1969), 
does  appear  to  be  slightly  more  sensitive  to  the  sea  surface  temperature 
variations.  The  similarity  of  the  results  is  closely  related  to  the  criti- 
cal relative  humidity  specification,  since  both  boundary  layer  versions 
have  large  evaporation  rates  which  necessitate  the  parameterization  of  con- 
vective  clouds.   In  this  sense  the  external  parameters  are  dominant  and 
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other  measures  must  be  used  to  compare  the  internal  parameters  of  the  two 
models. 

Model  Parameter  Distributions 

One  of  the  objectives  of  the  model  formulation  was  to  determine  equili- 
britun  atmospheric  conditions  for  given  ocean  temperatures,  rather  than 
simply  specifying  fixed  air-sea  temperature  and  moisture  differences.  Values 
of  sea-air  temperature  differences  are  shown  in  Figure  3,  for  the  CBL  with 
a  maximum  relative  humidity  of  955^.  Over  most  of  the  region  the  surface  air 
temperature  remains  essentially  constant.   A  somewhat  larger  sea-air  tempera- 
ture difference  is  predicted  in  more  intense  storms  that  occur  with  higher 
ocean  temperatures.  Near  the  central  region,  the  air  temperature  decreases 
rapidly.   In  the  minimum  hurricane  (27C)  the  temperature  decrease  is  to  a 
value  IC  below  the  water  temperature,  but  for  the  most  intense  (31C)  storm 
the  decrease  is  more  than  4C,  Riehl  (1963)  used  a  constant  difference  of 
2,5C  for  the  region  with  wind  speeds  exceeding  25  m  sec   ,   Bell  and  Tsui 
(1973)  found  the  mean  air  temperature  for  60  ship  observations  within 
100  n  mi  of  the  center  was  1,2C  less  than  the  ocean  temperature.  A  5C  dif- 
ference was  observed  as  the  center  of  hurricane  Hilda  passed  with  40  n  mi 
of  the  NOMAD  buoy  (see  Marcus  and  Smith,  1966).   Thus,  one  might  conclude 
that  the  air-sea  temperature  differences  from  the  model  are  not  inconsistent 
with  the  variation  determined  from  the  few  available  observations. 

A  number  of  other  parameters  were  examined  to  evaluate  the  differences 
between  the  model  results  with  the  two  boundary  layer  treatments.  Radial 
variations  of  several  of  these  parameters  are  illustrated  in  Figure  4  for 
two  cases  in  which  the  wind  speed  profiles  at  the  top  of  the  boundary  layer 
were  essentially  identical  (T  =  29C  in  Table  1). 
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In  the  CBL  model  the  stress  is  determined  from  the  iterative  solution 
for  the  friction  velocity  (u^) »  whereas  in  the  bulk  transport  equation  the 
stress  is  primarily  related  to  the  tangential  wind  speed  at  the  top  of  the 
boundary  layer.  Nevertheless,  the  surface  shearing  stress  (Tq)  profiles 
(Figure  4)  were  very  similar  for  both  models.   These  values  may  be  compared 
to  Tg  estimates  by  Miller  (1962)  from  data  collected  in  hurricane  Helene. 
Miller  calculated  the  shearing  stress  necessary  to  offset  the  inward  momen- 
tum transport  by  the  radial  wind  component,  which  was  estimated  from  air- 
craft data  at  800  mb  and  200  mb.   Average  T  values  were  60,  51  and  28  dynes 

-2 
cm   for  the  radial  increments  0-37,  37-74,  74-111  km,  respectively.  Two 

factors  may  explain  these  large  t  values  compared  to  those  from  the  model. 

Miller  assumed  the  surface  tangential  wind  component  was  equal  to  the  800  mb 

value,  while  a  reduction  of  0.7  was  used  in  the  model.   In  addition,  the 

radial  wind  component  estimates  by  Miller  were  larger  than  the  integrated 

V  values  from  the  model,  shown  in  Figure  4.  Thus  the  transfer  to  the 

ocean,  or  the  associated  radial  divergence  of  momentum,  was  smaller  for  the 

model  results.   Experiments  in  which  the  drag  coefficient  (and  thus  the 

-3 

heat  and  moisture  exchange  coefficients)  was  constant  (C  =  2.5  x  10  ) 

did  not  result  in  different  maximum  wind  speeds.   Even  though  the  heat  and 

moisture  fluxes  at  outer  radii  were  increased  in  these  cases,  the  6  at  the 

'      e 

eyewall  was  largely  unaffected.   Because  the  maximum  wind  speed  is  deter- 
mined by  the  empirical  relations  (2)  -  (4),  the  primary  effect  of  the  larger 
value  of  C  was  to  increase  the  frictional  stress  at  the  ocean  surface,  and 
this  was  balanced  by  increased  radial  advection  of  momentum.   In  this  respect, 
the  model  differs  from  others  (e.g.  Rosenthal  1971,  Ooyama  1969)  in  which 
the  ultimate  intensity  reached  by  the  model  storm  varies  inversely  with  C  . 

Another  indicator  of  the  storm  intensity  is  the  amount  of  precipitation. 
In  a  system  such  as  the  hurricane,  almost  all  of  the  moisture  in  the  inflow 
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layer  is  precipitated,  since  the  upper  level  outflow  contains  little 
moisture.   An  estimate  of  the  precipitation  is  then  the  vertical  moisture 
flux  through  the  top  of  the  boundary  layer  by  both  the  mean  and  convective 
motions.   The  moisture  source  due  to  both  radial  advection  and  evaporation 
is  included  in  this  precipitation  estimate.  An  integrated  radial  component 
of  3  m  sec   through  the  100  km  radius  with  a  moisture  content  of  20  gm  kg"''' 
contributes  about  12  cm  day   to  the  area-averaged  precipitation  within 
that  radius.   Since  evaporation  within  100  km  of  the  center  adds  only 
another  1  cm  day   to  the  precipitation,  the  total  is  closely  related  to 
the  strength  of  the  inflow.   It  might  be  noted  that  more  than  92^  of  the 
vertical  moisture  flux  was  by  the  mean  vertical  motion.  Maximum  precipita- 
tion rates,  which  occurred  at  the  eyewall,  for  the  bulk  transport  and  the 
CBL  models  compared  in  Figure  4  were  65  and  30  cm  day"  ,  respectively. 

Heat  Extraction  from  Ocean  -  " " 

Because  the  model  is  eventually  to  be  used  to  simulate  changes  in  ocean 
thermal  structure  changes  due  to  the  passage  of  a  tropical  cyclone,  the  total 
heat  extraction  and  its  radial  distribution  are  important  considerations. 
Black  and  Mallinger  (1972)  have  recently  summarized  various  estimates  of 
the  heat  flux  from  the  ocean.  These  estimates  were  based  on  bulk  transport 
equations  with  surface  temperatures,  moistures  and  winds  extrapolated  from 
aircraft  measurements.  Maximum  fluxes  occurred  near  the  radius  of  maximum 

winds,  with  the  estimates  ranging  from  less  than  800  to  more  than  6000  cal 

-2    -1 
cm   day  .  As  Black  and  Mallinger  point  out,  in  many  of  these  flux  calcu- 
lations the  sea  surface  temperatures  were  assvmied  equal  to  pre-storm  values, 
and  thus  may  overestimate  the  flux.  A  second  feature  of  these  calculations 
was  the  rapid  decrease  in  heat  flux  with  increasing  radius.  Nearly  all  of 
the  flux  estimates  at  five  times  the  radius  of  maximum  winds  were  less  than 
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one-half  of  the  peak  values.   Similar  variations  of  total  heat  flux  are 
evident  in  the  model  results  shown  in  Figure  4.  Maximum  values  occur 
within  a  few  km  of  the  radius  of  maximxim  winds,  but  are  considerably  smaller 
than  estimates  based  on  empirical  data.  The  values  in  Figure  4  represent 
the  total  heat  removed  in  the  first  24  hours  of  the  calculation,  with  the 

storm  stationary  and  the  ocean  temperature  fixed.  Maximum  values  of 

-2 
instantaneous  latent  and  sensible  heat  fluxes  at  24  hours  were  1.63  x  10 

-2       -2    -1 
and  0.31  x  10   cal  cm   sec   .  Thus  the  Bowen  ratio  near  the  eyewall 

was  about  20^.   By  contrast,  near  300  km  the  Bowen  ratio  was  approximately 

5^.   It  should  be  noted  that  the  heat  and  moisture  exchange  coefficients 

were  equal  to  2/3  C  in  the  model  with  the  bulk  transport  equations,  since 

it  was  expected  that  the  molecular  exchange  of  heat  and  moisture  at  the 

interface  would  be  inhibited  compared  to  the  momentum  exchange  (Robinson, 

1966).  A  series  of  experiments  was  also  rvin  with  the  heat  and  moisture 

coefficients  equal  to  the  C  value.  Latent  and  sensible  heat  fluxes  were, 

of  course,  increased  by  at  least  50^.   Extreme  values  of  heat  removal 

-2    -1 
ranged  from  650  to  6200  cal  cm   day   for  ocean  temperatures  of  27  and  31C, 

respectively.  However,  the  differences  in  the  temperature  fields  were 

less  than  0.5C,  with  warming  at  large  radii  and  enhanced  cooling  near 

the  center.  The  maximum  wind  speeds  were  increased  by  only  3-4  m  sec   over 

the  values  shown  in  Table  1.   Associated  increases  in  Tq,  v  and  maximiim 

y  r 

precipitation  rate  were  also  noted.   It  thus  appears  that  the  intensity  of 
the  storm  produced  by  the  bulk  transport  model  is  rather  insensitive  to  ex- 
change coefficients.   Specifying  coefficients  greater  than  the  minimum 
values  necessary  to  establish  the  66  possible  for  a  particular  ocean  tem- 
perature produces  proportionately  larger  heat  extraction,  but  not  signifi- 
cantly more  intense  storms.   In  future  experiments  with  the  ocean  surface 
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temperatures  predicted  by  the  heat  balance,  the  variation  of  C  might  be 
expected  to  have  a  larger  effect.   Perhaps  one  method  for  determining 
these  exchange  coefficients  will  be  to  determine  the  oceanic  heat  loss 
from  expendable  bathythermograph  traces  as  Black  and  Malliiiger  (1972)  have 
done  for  Hurricane  Ginger.   Calculation  of  the  oceanic  heat  budget  may  be 
more  reliable  than  the  use  of  bulk  transport  equations  with  assumed  ex- 
change coefficients  and  with  surface  winds,  temperatures  and  moistures 
extrapolated  from  aircraft  data. 

The  total  heat  flux  after  24  hours  in  the  CBL  model  was  generally 
smaller  than  that  calculated  in  the  bulk  transport  version  (see  Figure  4). 
Although  accurate  observations  of  the  heat  flux  within  the  hurricane 
boundary  layer  are  not  available,  the  heat  fluxes  from  the  two  models 
appear  to  be  underestimates,  especially  at  outer  radii.   A  lower  heat 
flux  may  have  been  a  result  of  the  assimiption  that  the  total  heat  flux 
was  proportional  to  the  density  (or  virtual  temperature) gradient  near 
the  air-water  interface.   In  the  mid-latitudes  the'  stability  is  primarily 
determined  by  the  large  air-sea  (or  air-land)  temperature  differences, 
and  the  moisture  content  is  generally  not  thought  to  be  important.  How- 
ever, in  the  tropics  significant  latent  heat  exchange  can  occur  with 
rather  small  air-sea  temperature  differences.  Assuming  similarity  between 
the  wind  and  virtual  temperature  profiles  allows  heat  transfer  even  if  the 
air  and  sea  temperatures  are  equal,  but  it  may  contribute  to  an  under- 
estimate of  the  total  heat  flux.  Another  factor  in  the  heat  flux  relation 
[see  (16)]  is  the  ratio  (a,  )  of  the  eddy  conductivity  (K^^)  to  the  eddy 
viscosity  (K^),  which  is  constant  if  the  profiles  are  similar.  However, 
the  value  of  the  ratio  over  the  sea  is  not  well  known.  Thus  it  has  been 
considered  a  disposable  parameter  in  the  model.   In  a  sense  this  is  analogous 
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to  the  experiments  described  above  in  which  the  heat  and  moisture  coeffi- 
cients were  varied  relative  to  C  .   Another  explanation  for  the  low  heat 
flux  estimates  could  be  that  the  u^  or  stress  values  are  underestimated. 
The  key  relation  for  the  surface  roughness  parameter  as  a  function  of  u^ 
may  also  be  varied  to  increase  the  frictional  drag  at  the  surface.  This 
was  the  second  approach  used  to  test  the  two-layer  model.   Point  values 
of  several  parameters  are  provided  in  Table  2  from  the  runs  with  T  =  29C. 
The  surface  roughness  coefficient  (C2  =  0.00428)  in  (13)  and  the  ratio  of 
heat  conductivity  to  eddy  viscosity  were  also  varied  in  runs  for  the  other 
sea  surface  temperatures,  and  the  resulting  changes  were  similar. 

Enhancement  of  the  heat  exchange  relative  to  momentum  exchange  (Exp.B) 
increased  the  maximum  wind  speed  about  4  m/sec.  This  increase  in  maximum 

wind  speed  ranged  from  3  to  7  m/sec  for  a  T  range  of  27  to  31C.  As  in 

w 

the  case  of  the  bulk  aerodynamic  model,  the  increase  in  surface  heat  ex- 
change beyond  a  certain  critical  value  does  not  dramatically  change  the 
predicted  intensity  of  the  storm  for  a  given  T  .  The  need  for  observational 
data  to  determine  the  proper  K^/K„  ratio  is  again  evident.   Increasing  the 
surface  roughness  in  the  model  (Exp.C)  also  resulted  in  more  intense  storms, 

with  the  increase  in  intensity  ranging  from  2  to  4  m/sec  for  T  values  be- 

w 

tween  27  and  31C.  As  indicated  earlier,  increases  in  the  surface  stress  in 
this  model  tends  to  increase  the  supply  of  energy  through  radial  advection 
without  greatly  retarding  the  wind  speeds  at  the  top  of  the  boundary  layer. 

A  general  increase  in  the  values  of  the  parameters  in  Table  2  may  be 
expected  with  the  increase  in  Vq,  for  example,  the  larger  radial  component. 
However,  in  Exp.  C,  v  was  increased  by  about  35^  over  the  control  value, 

compared  to  a  10'^  increase  in  Exp.  B.  The  stress  was  about  40^  larger  with 

-3 
an  implied  drag  coefficient  of  6.4  x  10   near  the  center  in  Exp.  C.   Even 
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at  300  km  the  drag  coefficient  was  2.0  x  lO"  .  The  associated  inflow 
angles  (ranging  from  14  to  27  degrees)  and  ratio  of  the  velocity  at  the 
top  of  the  surface  layer  to  that  at  the  top  of  the  boundary  layer  near 
the  center  were  also  significantly  changed  from  the  control  run  values. 
By  contrast,  the  above  factors  related  to  the  predicted  wind  fields  in 
Exp.  B  were  not  changed  beyond  the  expected  variations  due  to  the  in- 
creased maximum  wind  speed. 

The  total  heat  removed  from  the  ocean  in  Exps.  B  and  C  was  larger 
than  in  the  control  rxin.   As  expected,  the  increase  was  about  50^  when 

the  K^/Kj,  ratio  was  increased  from  2  to  3.   Peak  values  ranged  from  410 

-2    -1 
and  3940  cal  cm   day   for  T  equal  to  27  and  31C,  respectively.   En- 
hancement of  the  momentum  exchange  at  the  surface  produced  correspondingly 

larger  heat  fluxes,  with  the  peak  values  ranging  from  420  to  greater  than 

2    -1 
4670  cal  cm  day   for  the  same  T  values.   The  larger  precipitation  rate 

near  the  center  in  Exp.  C  is  evidently  due  more  to  a  larger  v  value 
than  to  the  increase  in  heat  gain,  because  the  50^  larger  heat  gain  in 
Exp.  B  is  accompanied  by  only  a  20^  increase  in  precipitation  rate.  This 
is  consistent  with  the  primary  moisture  source  being  at  larger  radii. 
The  inflow  angles  and  implied  drag  coefficients  in  Exp.  C  appear  excessive, 
particularly  at  outer  radii  where  more  observations  are  available  for  com- 
parison.  The  mechanism  through  which  the  heat  transfer  is  enhanced  is 
probably  in  the  wind-generated  spray  (Riehl,  1954),  which  exposes  a  much 
greater  surface  area  for  heat  transfer  from  the  water  to  the  air.   It 
appears  that  the  momentum  exchange  is  not  correspondingly  enhanced  by  the 
blowing  spray  (see  Monahan  1966,  and  Wu  1973).   If  this  is  the  case,  the 
ratio  of  K^/K»,  may  be  a  function  of  the  amount  of  spray,  or  the  surface 
roughness.  This  could  be  incorporated  in  the  model,  and  would  likely 
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increase  the  maximum  wind  speed  a  small  amount.   A  different  representa- 
tion of  the  heat  loss  from  the  ocean  versus  radius  would  also  result. 
The  flexibility  of  the  CBL  model,  and  the  fact  that  future  research 
results  from  air-sea  interaction  studies  will  likely  be  cast  in  the 
framework  of  similarity  theory,  or  external  parameters  such  as  Z  ,  is 
thought  to  be  sufficient  to  justify  this  more  complex  and  time-consuming 
boundary  layer  model. 

CONCLUDING  REMARKS 

Summary 

An  empirical  wind  profile  and  an  expression  for  the  maximum  wind 
speed  in  terms  of  the  radial  increase  in  0  have  been  used  to  specify  the 
winds  at  the  top  of  the  boundary  layer  in  a  sjnnimetrical  hurricane.  Using 
these  winds  with  the  bulk  aerodynamic  transport  equations,  and  as  an  input 
to  a  two-layer,  baroclinic  marine  boundary  layer  model  by  Cardone  (1969), 
solutions  have  been  obtained  for  the  moisture  and  potential  temperature 
distributions  within  the  boundary  layer.  The  vortex  in  turn  responds  to 
the  equivalent  potential  temperature  gradient  induced  by  the  balance  be- 
tween the  heat /moisture  sources  and  the  advection  by  the  boundary  layer 
winds.   Equilibrium  maximum  wind  speeds  vary  from  minimal  hurricane  force 
to  80  m  sec   for  ocean  temperatures  between  27C  and  31C.  This  might  be 
considered  the  maximum  storm  for  the  specified  ocean  conditions;  in  actual 
storms  the  external  environmental  effects  would  determine  if  these  wind 
speeds  would  be  achieved.   For  a  given  ocean  temperature,  maximum  wind 
speed  in  the  model  is  shown  to  be  very  sensitive  to  the  maximum  relative 
humidity  of  the  surface  air  near  the  eyewall,  which  is  specified  as  an 
external  parameter  in  an  attempt  to  parameterize  the  convective  flux  through 
the  top  of  the  boundary  layer.  The  bulk  aerodynamic  and  the  Cardone 
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boundary  layer  approach  give  similar  maximum  wind  speeds  for  the  same 
external  conditions.   There  is  a  certain  arbitrariness  in  the  total  heat 
flux  from  the  ocean  for  a  given  viscosity  in  both  boundary  layer  versions. 
Observations  of  heat  loss  from  the  ocean  are  needed  to  establish  the 
appropriate  constants  in  the  high  wind  speed  conditions  that  characterize 
hurricanes.   Various  parameters  predicted  by  the  Cardone  model,  such  as 
the  radial  mass  flow  and  inflow  angles,  depth  of  the  boundary  layer,  verti- 
cal motion  and  turbulent  fluxes,  may  be  evaluated  by  aircraft  equipped 
with  inertial  guidance  systems  that  are  available.  -  _..'.-■] 

Possible  Applications  of  the  Model 

The  original  objective  was  to  develop  an  atmospheric  vortex  capable 
of  simulating  the  effect  of  a  hurricane  on  the  ocean  thermal  structure. 
Since  the  ocean  temperature  affects  the  intensity  of  the  vortex,  a  mutual 
interaction  may  occur  if  significant  cooling  occurs  in  the  ocean.   Leipper 
and  Volgenau  (1970)  suggested  that  the  subsurface  thermal  structure  must 
also  be  considered  in  forecasting  intensity  changes  in  hurricanes.  Many 
experiments  are  planned  to  evaluate  the  effect  of  gradients  in  both  sea 
surface  temperature  and  thermocline  depth.  A  layer  ocean  model  similar  to 
that  used  by  Kraus  and  Turner  (1967)  is  to  be  tested. 


The  adjustment  time  in  the  model  is  relatively  short  (see  Figure  5). 

For  T  =  31C,  the  increase  in  Vq    from  30  to  70  m  sec   occurred  within 

w  Umax 

11  hrs.   Thus  the  model  experiments  can  be  run  with  a  variety  of  sea  tem- 


perature conditions.   In  the  present  form  the  vortex  is  stationary,  so 
that  allowing  the  sea  temperature  to  vary  will  simulate  the  rate  of  decay 
of  the  vortex  due  to  cooling  at  inner  radii.   The  method  of  moving  the 
vortex  will  be  based  on  the  integrating  effect  of  strong  tangential  winds 
in  producing  properties  which  are  nearly  axisymmetric.   Application  of  the 
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two-dimensional  model  along  four  or  more  diagonals  will  be  used  to  esti- 
mate the  symmetric  component --which  determines  the  intensity  of  the  storm. 
In  this  way  the  model  vortex  may  be  moved  over  a  varying  ocean  thermal 
field.   One  advantage  of  the  present  model  over  previous  ntimerical  models 
is  the  possibility  of  simulating  the  symmetric  component  of  actual  storms. 
The  model  can  be  initialized  to  a  particular  maximum  wind  speed  by  adjust- 
ing the  exchange  coefficients  or  other  external  parameters,  for  the 
specified  ocean  temperature.   Experiments  are  planned  to  test  this  capa- 
bility with  real  data  by  moving  the  model  storm  along  the  actual  track  with 
the  external  parameters  held  constant.   One  should  then  be  able  to  assess 
the  predicted  changes  in  maximum  intensity  of  the  storm  and  in  the  ocean 
thermal  structure. 

Finally,  the  model  could  be  useful  in  testing  proposed  hurricane  modi- 
fication schemes  that  are  dependent  on  reduction  of  the  interface  fluxes 
of  heat  and  moisture,  or  the  inflow  within  the  boundary  layer.   It  has 
been  proposed  (Project  Stormfury  1968)  that  spreading  of  mononuclear  films 
could  inhibit  evaporation  from  the  ocean  and  thus  reduce  the  energy  source 
for  hurricanes.   Results  from  the  model  verify  this  reasoning,  particularly 
if  the  moisture  content  reduction  could  be  maintained  in  the  inner  region 
of  the  storm.   With  currently  available  organic  materials,  the  film  is 
dispersed  by  the  high  wind  speeds.   It  might  be  possible  to  efficiently 
test  many  patterns  of  material  spreading  to  determine  the  most  effective. 
There  may  be  some  optimum  radius  beyond  which  the  materials  cannot  be 
economically  laid  down  for  a  given  benefit  from  reducing  the  maximum  wind 
speed.   A  more  recent  proposal  by  Gray  (1973)  to  modify  the  hurricane  by 
carbon  dust  seeding  within  the  boundary  layer  might  also  be  tested  with 
this  model.   Gray's  hypothesis  is  that  artificial  warming  of  the  atmosphere 
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boundary  layer  due  to  absorption  of  insolation  by  the  carbon  dust  may  cause 
an  additional  mass  flux  out  of  the  boundary  layer.   A  decrease  in  the  low- 
level  inflow  is  expected  to  decrease  the  inner  core  maximum  wind  velocities, 
and  thus  reduce  the  storm  damage.   The  potential  benefit  of  the  various 
hurricane  modification  schemes  is  enormous— and  so  also  is  the  possible 
damage  from  misguided  attempts.   Wherever  possible,  these  schemes  should 
be  tested  on  a  variety  of  numerical  models,  including  the  type  proposed  here. 
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TABLE  1.   Maximiun  wind  speed  (m^ec)  in  the  model  as  a  function  of  sea 

temperature  for  the  two  boundary  layer  versions.  The  critical 
moisture  distribution  was  by  specification  of  the  cloud  base 
or  by  the  maximum  relative  humidity.   Values  in  parentheses 
indicate  speed  if  environmental  0e  values  exceeding  350K  are 
used  in  E.q.  (5)  rather  than  a  fixed  6e  =  350K  at  r  . 

Ocean            Bulk  Aerodynamic  CBL 

Temperature      Moisture  Specification  Moisture  Specification 

(deg  C) Cloud  base     75  to  95^ Cloud  Rase     75  to  95^ 

31  70  76  65(51)  76(62) 

30  58  64  56(47)  66(61) 

29  48  56  44(43)  56 

28  36  47  32  46 

27  36  33 
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TABLE  2.   Significant  parameters  at  inner  (I),  middle  (M)  and  outer  (0) 
regions  from  CBL  model  upon  variation  of  the  surface  roughness 
(coefficient  C2)  and  the  ratio  of  heat  conductivity  to  eddy 
viscosity  (a,  ) •   Ocean  temperature  held  fixed  at  29C. 


Variable 
C2 

^  -1 

Vnfm  sec   ] 

V  [m  sec   J 
r 


-2- 


Stress  [dyne  cm  ] 
V,. 


Inflow 
Angle  [deg] 
Cjj  X  10-^   •;, 

Precip  [cm  day  ] 
Total  heat  removed 
[cal  cm   day  ] 


Region 


I 

I 
I 
I 
0 
I 
0 

I 

0 

I 
I 

M 


A 

Experiment 
B 

C 

0.00428 

0.00428 

0.025 

2 

3 

2 

56.5 

60.5 

58 

1.4 

1.5 

1.9 

37 

43 

50 

.60 

.62 

.50 

.85 

.83 

.82 

22 

23 

27 

10 

10 

14 

2.9 

3.1 

6.4 

1.1 

1.0 

2.0 

33 

38 

51 

980 

1460 

1670 

180 

100 

115 
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LIST  OF  ILLUSTRATIONS 

Figure  1.    Schematic  model  of  atmospheric  vortex  and  associated  Inflow 

layer.  The  exponent  x  is  -1  v/ithin  r.  and  1/2  between  r. 

and  r   (see  text) . 
o 

Figure  2.    Schematic  of  the  Cardone  (1969)  boundary  layer  model. 

Figure  3.    Radial  variation  of  the  air  minus  sea  temperatures  for  the 

CBL  with  the  ocean  temperature  fixed  at  27,  29  and  31C.  The 
curves  indicate  the  steady  state  values  with  the  ocean  tem- 
perature held  fixed. 

Figure  4.    Comparison  of  radial  profiles  of  selected  parameters  from 
the  bulk  aerodjmamic  boundary  layer  equations  (solid  line) 
and  the  Cardone  (CBL)  model  (dashed).  The  radial  velocity 
at  the  top  of  the  surface  layer  (v*) ,  from  the  CBL  model  is 
shown  in  the  middle  panel. 

Figure  5.    Variation  of  maximum  wind  speed,  starting  from  an  initial 

wind  speed  of  30  m  sec   ,  for  different  ocean  temperatures. 
The  arrows  indicate  the  time  at  which  the  specific  hiamidity 
value  at  the  eyewall  first  exceeded  the  critical  relative 
humidity. 
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Figure  1.    Schematic  model  of  atmospheric  vortex  and  associated  inflow 
layer.   The  exponent  x  is  -1  within  r.  and  1/2  between  r- 
and  r   (see  text) . 
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Figure  3,    Radial  variation  of  the  air  minus  sea  temperatures  for  the 

CBL  with  the  ocean  temperature  fixed  at  27,  29  and  31C.   The 
curves  indicate  the  steady  state  values  with  the  ocean  tem- 
perature held  fixed. 
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Figure  4.    Comparison  of  radial  profiles  of  selected  parameters  from  the 
bulk  aerodynamic  boundary  layer  equations  (solid  line)  and 
the  Cardone  (CBL)  model  (dashed).   The  radial  velocity  at  the 
top  of  the  surface  layer  (v*),  from  the  CBL  model  is  shown  in 
the  middle  panel. 
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Variation  of  maximum  wind  speed,  starting  from  an  initial 
wind  speed  of  30  m  sec"^,  for  different  ocean  temperatures. 
The  arrows  indicate  the  time  at  which  the  specific  humidity 
value  at  the  eyewall  first  exceeded  the  critical  relative 
humidity. 
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